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The renormalized Dyson-Schwinger equation for the quark propagator is studied, in Landau gauge, 
in a novel truncation which preserves multiplicative renormalizability. The renormalization con- 
stants are formally eliminated from the integral equations, and the running coupling explicitly 
enters the kernels of the new equations. To construct a truncation which preserves multiplicative 
renormalizability, and reproduces the correct leading order perturbative behavior, non-trivial can- 
cellations involving the full quark-gluon vertex are assumed in the quark self-energy loop. A model 
for the running coupling is introduced, with infrared fixed point in agreement with previous Dyson- 
Schwinger studies of the gauge sector, and with correct logarithmic tail. Dynamical chiral symmetry 
breaking is investigated, and the generated quark mass is of the order of the extension of the in- 
frared plateau of the coupling, and about three times larger than in the Abelian approximation, 
which violates multiplicative renormalizability. The generated scale is of the right size for hadronic 
phenomenology, without requiring an infrared enhancement of the running coupling. 

PACS numbers: 12.38.Aw, 14.65.-q, ll.15.Tk, ll.10.Gh 



I. INTRODUCTION 



In the Standard Model of the strong, weak and electromagnetic forces, the interactions are quantitatively described 
by gauge held theories. Quantum chromodynamics is a non- Abelian gauge theory, and the proof of its renormaliz- 
ability M and discovery of ultraviolet asymptotic freedom || have been milestones in its acceptance as theory of the 
strong interaction. For large momenta the coupling becomes small, and perturbation theory seems an appropriate 
calculational tool. However, for small momenta the coupling grows large and adequate methods have to be used to 
study nonperturbative phenomena like confinement, dynamical chiral symmetry breaking and bound state formation. 
One such method is the study of the Dyson-Schwinger equations (DSE) |}| , and their phenomenological applications 
to hadronic physics is a subject of growing interest |^|. 

Dynamical chiral symmetry breaking can be studied by means of the quark propagator DSE, also called gap 
equation. The quark equation is part of an infinite tower of integral equations relating all the Green's functions of 
the quantum held theory, and a truncation is necessary to be able to solve it. The often used Abelian approximation 
introduces an effective running coupling in the kernels of the quark equation, based on features of QED. A number 
of studies of the gap equation were performed in this approximation with effective strong couplings that are infrared 
enhanced []5|, |^, Q] , infrared vanishing ^) , infrared vanishing with large enhancements in the intermediate region jl0| , 
and infrared hnite [ pT| [l^ ]. These studies all find that dynamical chiral symmetry breaking can be triggered, provided 
some parameter in the model exceeds a critical value, and a large integration strength is needed in the self-energy 
kernels to achieve large enough "constituent" quark masses, needed for hadron phenomenology The Ansatz of 
a strong infrared enhancement in the kernels of the quark DSE is therefore a major ingredient in phenomenological 
studies, and the search for the source of such an enhancement, coming from the gluon propagator, quark-gluon vertex 
or their combination is an important held of investigation. 

The kernel of the gap equation not only reflects the strong interaction between quarks and gluons, but also embodies 
the complicated structure of the QCD vacuum, which in the continuum contains self-interacting gluon helds as well 
as ghost helds. Understanding the infrared behavior of the gluon and ghost propagators, and of the running coupling 
is therefore of great importance for the study of dynamical mass generation |l3| . 

Early studies of the Dyson-Schwinger equation for the gluon propagator in the Landau gauge seemed to indicate that 
the gluon propagator could be highly singular in the infrared [14], [l5, 16, 1^ ], possibly providing the above mentioned 
infrared enhancement in the kernels of the gap equation. However, these gluon propagator studies neglected any 
contribution of the ghost fields, and required the ad hoc cancellations of certain leading terms in the equations. It is 
therefore far from certain that these solutions reflect the correct QCD infrared behavior. 

More recently, studies of the coupled set of Dyson-Schwinger equations for the gluon and ghost propagators |T^, 
E| have shown that the coupling of ghost and gluon fields plays a crucial role in the generation of a consistent 
infrared behavior of QCD. In the infrared the gluon and ghost propagators are power behaved, and the strong running 
coupling has an infrared fixed point. Moreover, the leading infrared power in the gluon vacuum polarization depends 
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only on the ghost loop, and not on the gluon loop nor any other diagram in the gluon equation. Although the precise 
values of the infrared powers depend on the vertex Ansatze and other details of the truncation, all these studies 
showed that the full ghost propagator is more divergent than its bare counterpart, while the full gluon propagator 
is less divergent. Recent lattice calculations also support the infrared propagator power laws and fixed point of the 
running coupling pl| , p2| , 23 . As neither the gluon propagator, nor the running coupling are infrared enhanced in 



the Landau gauge, one seems to be at loss to explain how enough strength can be achieved when breaking the chiral 
symmetry in the gap equation. The problem to reconcile the DSE results of the gauge sector with the need for a 
strong infrared enhancement in the effective coupling might hint to a deficiency in the Abelian approximation to 
the gap equation. A first attempt to extend the study of dynamical chiral symmetry breaking beyond the Abelian 
approximation, by solving the coupled set of DSEs for the quark, gluon and ghost propagators simultaneously 
underlines this conundrum by finding a generated fermion scale that is a factor two too small for phenomenological 
purposes. One can then wonder if the approximations, introduced to solve the gap equation, violate some essential 
property of quantum chromodynamics, such that the strength of dynamical chiral symmetry breaking is significantly 
weakened? And if so, can we construct a truncation which preserves this property, and would naturally generate a 
large enough dynamical quark mass? 

Multiplicative renormalizability (MR) is an important feature of gauge field theories, proven for perturbative renor- 
malization Jl[], but believed to hold for any renormalization scale, and it is often ruined by the approximations to the 
vertices and vertex renormalization constants introduced in DSE studies to truncate the infinite tower of equations. 
None of the truncations introduced in the DSE studies mentioned above, on the quark, gluon and ghost DSEs, preserve 
multiplicative renormalizability. 

In a recent paper |25|] we have discussed how the DSE studies [ffM [l9|, |2Cj of the gauge sector violate the multiplicative 
renormalizability of the gluon and ghost propagators, and we developed a new truncation scheme which preserves 
multiplicative renormalizability, and showed how this induced important changes in the equations. We investigated 
the consequences for the analysis of the infrared behavior of the propagators, and found that the ghost loop, gluon 
loop, 3-gluon and 4-gluon diagrams all potentially contribute to the leading infrared power in the gluon propagator 
DSE, which is at variance with the conclusions of Refs. Jl8|, [h| ^o|, where only the ghost loop contributed to leading 
order. Although the combination of ghost and gluon loops alone does not seem to allow infrared power solutions for 
the propagators [^5| , there are strong indications that the infrared contribution of the 4-gluon diagram, also called 
squint diagram, is such that infrared power solutions do exist when all the diagrams are taken into account ]2~q |. 

In this paper we construct an MR preserving truncation to the quark DSE using an approach similar to that 
presented in Ref. ^{| for the gauge sector, and investigate dynamical chiral symmetry breaking. We show that the 
preservation of multiplicative renormalizability in the truncated quark equation can naturally generate a phenomeno- 
logically acceptable mass scale, while being consistent with the properties of the gluon and ghost propagators known 
from the DSE studies of the gauge sector, without infrared enhancement of the strong running coupling. 

We first reformulate the quark equations, such that the Landau gauge equations are free of renormalization con- 
stants. In this formulation the self-energy kernels are explicitly driven by the running coupling, and the multiplicative 
renormalizability of the quark dressing function, and renormalization point invariance of the dynamical mass are 
manifest. This allows us to design a novel truncation scheme for these formal equations, which respects both the 
multiplicative renormalizability and the resummed perturbative limit of the solutions. The full vertex is genuinely 
nonperturbative in our truncation, and non-trivial cancellations in the self-energy integrals are assumed to strip the 
vertex from its complete dressing, leaving only bare vertices in the truncated equations. 

In the new truncated quark equations, the running coupling is the only object carrying information about the 
nonperturbative effects in the gauge sector. The back reaction of the quarks on the running coupling could in 
principle be included through a self-consistent study of the coupled quark-gluon-ghost DSEs. However, to perform a 
first quantitative calculation of dynamical chiral symmetry breaking, we introduce a model running coupling based 
on prior knowledge about its infrared fixed point, acquired from the ghost-gluon DSE p0| , p5| , the resummed leading 
order perturbative behavior, and the qualitative features of the transition between both regions . 

Dynamical chiral symmetry breaking is then investigated using this model running coupling in the MR truncation of 
the quark DSE. Although the equations in our new truncation look very similar to those in the Abelian approximation, 
the slight modifications ensuring MR yield important quantitative differences: the dynamically generated mass in the 
MR truncation is boosted by about a factor three compared to the Abelian approximation with same running coupling. 
Our new findings show that no infrared enhancement, nor fine-tuning of the strong running coupling is needed: a 
running coupling having a reasonably valued infrared fixed point (ao ~ 2.6), which smoothly bends over into the 
perturbative logarithmic tail, is able to generate masses allowing for realistic hadron phenomenology. 
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II. THE QUARK EQUATION 



Starting from the spinor quark DSE p], we derive two scalar equations for the mass function and renormalized quark 
dressing function, using the definitions of renormalized fields and coupling. The ensuing quark equations are part 
of an infinite tower of integral equations relating all the Green's functions of QCD, and a truncation is necessary to 
be able to solve them. We briefly discuss the Abelian approximation, which violates multiplicative renormalizability, 
before constructing a novel truncation which preserves the multiplicative renormalizability of the quark solutions. 

The Dyson-Schwinger equation for the quark propagator in QCD, formulated in Euclidean space, is 



A 

(27T) 



IT' (p, q,r)S F (q) (q, p, -r) (r) , 



4 M 



(1) 



where Sp and S F are the full and bare quark propagators, D^ v is the full gluon propagator, TJf' and T% 9 are the 
bare and full quark-gluon vertices, the color factor Cf — (AT? — l)/2N c — 4/3 for N c = 3, go is the bare coupling, 
and r = p — q. 

The most general expression for the full quark propagator can be written as 



Sf{p) 



Z(P 2 



ip-j + M(p 2 ) ' 



(2) 



where Z(p 2 ) is the quark dressing function, and M(jp ) is the mass function. The bare propagator for a quark with 
bare mass m is 



S F (P) = - 



1 



ip ■ 7 + mo 

The full gluon propagator in a general covariant gauge is given by 



Pp.Pv \ F{p 2 ) Pjjh^ 

..0 ' S _A 



(3) 



(4) 



where F{p 2 ) is the gluon dressing function, and for future use we also give the general form of the full ghost propagator: 



A I V G( P 2 ) 

a (p) = — zr~ 



(•5) 



where G(p 2 ) is the ghost dressing function. 

After substituting the propagators (0) - (^) in Eq. (|l|), we can derive a set of two coupled equations for the quark 
dressing function and the mass function by multiplying Eq. (Q) respectively by — ip ■ 7 and I4, and taking the trace in 
spinor space: 
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Z(p 2 ,A 2 ) 

M(p 2 ) 
Z(p 2 ,A 2 ) 

where 

U z (p 2 ,q 2 ,r 2 ,A 2 ) 
U M (p 2 ,q 2 ,r 2 ,A 2 ) 



= 1 + 
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C F 
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4r 2 



rTr 



q ■ 7 

M(q 2 ) 



T™{q,p-r,A 2 ) 



F{r 2 ,A 2 ) r 2 



F(r 2 ,A 2 



(6) 
(7) 

,(8) 
(9) 



with 5t±{f) — g^ v — r^r v /r 2 , and the kernels Uz and Um also depend implicitly on Z and M through the full quark- 
gluon vertex T1 9 . The vector and scalar quark self-energy integrals are herein regulated with an ultraviolet cutoff A, 
yielding a A 2 dependence of the regularized dressing functions, and the A 2 dependence of the bare mass TOo(A 2 ) is 
such that M(p 2 ) is finite and independent of A 2 . 
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Renormalization 

The full, regularized Green's functions are potentially divergent when the ultraviolet cutoff A is taken to infinity, 
and they are renormalized by applying the principles of multiplicative renormalization. The quark field is renormalized 
by 

Z(p 2 ,A 2 ) = Z 2 (» 2 ,A 2 )Z R (p 2 ,[x 2 ), (10) 

where Z R is the renormalized quark dressing function, Z2 is the quark field renormalization constant, and at the 
renormalization point Z R (fi 2 , fi 2 ) = 1. 

When addressing mass renormalization it is important to observe that each choice of bare mass parameter mo (A 2 ) 
defines a different physical theory p7| , and for QCD each quark flavor corresponds to a different value mg(A 2 ). The 
running mass function M(p 2 ) is renormalization point invariant as it is unambiguously determined by m (A 2 ). The 
bare mass parameter is chosen such that the running mass function M(jp ) is finite, and takes the value M(p 2 ) — m^, 
where the renormalized mass parameter is determined, in practice, by matching the calculated value of some mass 
dependent observable to its measured value. Note that we speak of mass parameter instead of mass, to emphasize 
that the quarks are confined and have therefore no pole mass. 

The renormalization of the quark-gluon vertex introduces a vertex renormalization constant Z\f, which is related 
to the renormalization of the strong coupling by 

Z lf (n 2 ,A 2 ) 

where g(p 2 ) is the value of the renormalized coupling at the renormalization scale. Gauge invariance of the renormal- 
ized theory ensures the universality of the renormalized coupling for the quark-gluon, ghost-gluon and triple-gluon 
interaction, such that also 

Z^V.A^G^A 2 ) 2 ^ 3 3/ V,A 2 ) 

9{p } = 5o(A } = w,a 2 ) 5o(A } ' (12) 

where Z\ and Z\ are the ghost-gluon and triple-gluon vertex renormalization constants. Z% and Z% are the renormal- 
ization constants for the gluon and ghost fields, 

F(p 2 ,A 2 ) = Z 3 (f, 2 ,A 2 )F R (p 2 ,^) , G( P 2 ,A 2 ) = Z 3 (n 2 ,A 2 )G R (p 2 ,[i 2 ), (13) 

where the renormalized gluon and ghost dressing functions are defined such that F r (/j?,ii 2 ) = G_r(/i 2 , /1 2 ) = 1. 
Now that all the necessary renormalized quantities have been defined, we multiply Eqs. (^), (0) by Z2 , such that 



their left hand sides become finite, and introduce the renormalized dressing functions using Eqs. (|10| ), (13), and the 
renormalized coupling g(^ 2 ) using Eq. (|ll|): 



Zr(p 2 ,H 



2\ 



Z 2 (^A 2 ) + ^4 2 )^(^,A 2 ) / A Z R (q 2 , fi 2 ) g ff ^ } U z {p 2 , q 2 , r 2 , A 2 ) , (14) 

Z^pf^ = Z 2 (n 2 ,A 2 )m (A 2 ) -^ a ^ 2 )Z 2 f ^ 2 ,A 2 ) J d\ Z R (q 2 , M 2 ) M(q 2 ) g ff ^ } U M (p 2 , q 2 ,r 2 , A 2 ) , 

(15) 

where we defined a(/j, 2 ) = g 2 (/i 2 ) / Air . Note that the integrals in the renormalized equations are preceded by a 
renormalization factor Z 2 ^(/i 2 , A 2 ), originating from the renormalization of the coupling. Although it is customary 
to absorb one Z\ factor in the full unrenormalized vertex in order to renormalize it, we deliberately do not proceed 
like this here, so that the full vertex depends on the momenta flowing in the vertex, and on A, but not on \i. The 
self-energy integrals have a A dependence consistent with perturbation theory, which is cancelled by the seeds of Eqs. 
(|l4|), ( fl5|) to yield a finite dressing function Z R and mass function M. 

Abelian approximation 

A truncation often used in previous studies to decouple the quark equations (pT[), ( |i~5| ) from the infinite tower of 
DSEs is the bare vertex, Abelian approximation. There, ghost field contributions are neglected, and one introduces 
the substitution 

«( M 2 ) z 2 f F R (r 2 , /x 2 ) r«"° r« — > « C V) r<j?>° r™>° , (ie) 



5 



based on QED gauge invariance properties, and on the assumption that the bare vertex approximation might be 
appropriate in the Landau gauge. In this approximation a cS (r 2 ) can be considered as an effective running coupling. 
In the Abelian approximation the quark equations ( 53 ) , (|TJ|) become 



Z^) - *<^AVg / ^W-^^M, (17) 
M (P 2 ) A 2x_ C F f^_„ a off (r 2 ) 



Z 2 (^,A 2 ) TO0 (A 2 ) -%\d\ Z R (q 2 ^ 2 ) Mtf) ^^ UbtfrfS), (18) 



where t/^ and are calculated by replacing the full, regularized, unrenormalized quark-gluon vertex T 99 by the 
bare vertex, T 99 ' = i<y v , in Eqs. (|), (§). 

It is easy to convince oneself that Eqs. (|l7|), (|l8| ) of the Abelian approximation violate multiplicative renormaliz- 
ability, and this will be illustrated further in the result section. 



MR truncation 



The aim of this paper is to derive a novel truncation scheme, which preserves the multiplicative renormalizability of 
the solutions, starting from the exact renormalized equations (fl4|), jl5|). This will be achieved by specific manipulations 
of the renormalization constant Z\f, and by assuming non-trivial cancellations involving the full quark-gluon vertex 
in the quark self-energy loop, as specified in the remaining of this section. 

The nonperturbative expression for the quark-gluon vertex renormalization constant Z\ / is not known, but we can 
eliminate it from our equations using Eqs. (|Tl|), (|l2|), which state the universality of the strong coupling: 



Z lf (fi 2 ,A 2 



W,a 2 ; 



^( M 2 ,A 2 ). 



(19) 



It is interesting to note the similarity between the QCD relation (|19|) and the QED analog, Zif = Z 2 , in a theory 
without ghost fields. 

As Z\ f is independent of the integration momentum, the factor Z 2 ^ can be moved inside the self-energy integrals 
of Eqs. ( |T4| ) , (p^| ), and factorized in a A and a \jl dependent part by applying Eq. (|l9|), and consequently using the 
definitions (|l0|), (|l3|) to replace Z2 and Z3 by ratios of unrenormalized to renormalized dressing functions at momenta 
of our choice: 



7 A ^ Z{g 2 ,A 2 )G R {r 2 ^ 2 ) ~ ? (u 2 ,2, 



(20) 



where we have chosen to use the quark momentum of the self-energy loop to rewrite Z2, and the gluon momentum 
for Z 3 . After substituting the square of Eq. ( gtj ) in the quark equations (|l4|), ( p"5| ) we find 
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M(p 2 ) 
Zr(p 2 ^ 2 ) 



Z 2 ( A1 2 ,A 2 ) + ga( A1 2 )Z 1 2 ( M 2 ,A 2 ) 



1 



F R (r 2 ^ 2 )G 2 R (r 2 ,„ 2 ) 



Z R {q\^) g 2 + M 2 ( g 2 ) 



Z 2 ( A/ 2 ,A 2 )m (A 2 



Cf_ 
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a{^)Z^,A 2 ) 



M(q 2 ) F*(r 2 , M 2 )G 2 (r 2 , M 2 ) 



Zn(q 2 ^ 2 ) 



M 2 {q 2 ) 



Z 2 (q 2 ,A 2 ) 
G 2 (r 2 ,A 2 ) 

^V,A 2 ) 
G 2 (r 2 ,A 2 ) 



U z (p 2 ,q 2 ,r 2 ,A 2 ) 



U M ( P \q 2 ,r 2 ,A 2 ) 



From Eqs. (jl2|), ( |l3| ) it is easy to see that the product 

a(q 2 ^ 2 ,A 2 )^a^ 2 )Z 2 1 ( i , 2 ,A 2 )F R (q 2 ,^)G 2 R (q 2 ,^ 2 ) 
can also be written in terms of unrenormalized quantities only: 

a(q 2 ,v 2 ,A 2 ) = 9 -^f±F(q 2 ,A 2 )G 2 (q 2 ,A 2 ), 

47T 



(21) 



(22) 



(23) 



(24) 
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which shows that a(q 2 , /i 2 , A 2 ) is renormalization group invariant, i.e. independent of /j, 2 . Therefore, a(q 2 , /i 2 , A 2 ) = 
a(q 2 , q 2 , A 2 ), and taking into account the renormalization conditions for Fr and G R , and the fact that Z\ = 1 in the 
Landau gauge f28| , we see that, in that gauge, the product a(q 2 ,fi 2 ,A 2 ) defined in Eq. (23) is nothing else but the 
running coupling a(q 2 ) (see also Refs. JT4], |l^, ff9|). One can identify the product (^3|) in Eqs. (pl|), (p2|), and in the 
Landau gauge these equations become: 



1 



Zr{j> 2 ^ 2 ) 

M{p 2 ) 
Zr(p 2 ^ 2 ) 



= Z 2 ( M 2 ,A 2 )m (A 2 )-g Jd\ 



a(r 2 ) 



Z R {q 2 ^ 2 ) q 2 + M 2 (q 2 ) 



Z 2 (q 2 ,A 2 



G 2 (r 2 ,A 2 ) 



U z (p 2 ,q 2 ,r 2 ,A 2 ) 



a(r 2 ) 



M{g 2 ) 

Z R {q 2 ^ 2 ) q 2 + M 2 (q 2 ) 



Z 2 (q 2 ,A 2 ) 



G 2 (r 2 ,A 2 



U M (p 2 ,q 2 ,r 2 ,A 2 ) 



(25) 
(26) 



In numerical calculations the seeds of the equations will be eliminated in the usual way by subtracting each equation 
at two different momenta, and imposing the renormalization conditions. Note that, except for the seeds of the 
integral equations (p5|), (P6|), the renormalization point dependence only enters through the renormalized quark 
dressing function Z R , so that the multiplicative renormalizability of the dressing function, and renormalization point 
invariance of the mass function are manifested in a very clear way by the 1/Z R structure on both sides of the equations. 
Furthermore these equations are still exact as no approximations have been introduced in their derivation, and we 
have shown how the running coupling explicitly enters the self-energy kernels, and drives the mass function, in a 
Non-Abelian gauge theory. 

The kernels Uz and Um, defined in Eqs. (^), (JoJ) , contain the full, unrenormalized quark-gluon vertex, thus coupling 
the propagator equations (^5|), (|26| ) to the higher order DSEs, and in order to perform a numerical calculation we 
introduce an approximation which decouples the quark equations from the vertex DSE. 

A simple truncation which preserves the MR properties of the solutions, and reproduces the leading order renor- 
malization group equation (RGE) improved perturbative behavior consists in assuming that the factors Z 2 /G 2 in 
the square brackets of Eqs. (|2^), (|2^) cancel both the nonperturbative vertex and loop corrections in the self-energy 
integrals. This truncation is analogous to the one introduced in the DSE for the gluon propagator in Ref. j2^]. The 
kernels between square brackets in Eqs. (|25|), (E6J) are thus approximated by 



Z\q 2 .A 2 \ 

G 2 (r 2 ,A 2 ; 
Z 2 (g 2 ,A 2 ) 
G 2 (r 2 ,A 2 ) 



■U z (p 2 ,q 2 y,A 2 ) 
U M (p 2 ,q 2 ,r 2 ,A 2 ) 



(27) 
(28) 



where 17% and are calculated by replacing the full, regularized, unrenormalized quark-gluon vertex T* 9 by the 
bare vertex, r« 9 <° = i<y v , in Eqs. (§), (§). After substituting the Ansatze @, (§|) in Eqs. ©, @ we find 
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Zr(p 2 ^ 2 ) 
M(p 2 ) 



^A 2 ) + g/d 4 



a(r 2 ) 



Z R {q 2 ^ 2 ) q 2 + M 2 (q 2 ) 



Z R (p 2 ^ 2 ) 
where the Landau gauge kernels (£ 



Z 2 ( Ai 2 ,A 2 )m (A 2 ) 



Cf_ 

47T 3 



0) are given by 
1 



U z (p 2 ,q 2 ,r 2 ) 



d*q 



3p-q 



a(r 2 ) 



M{g 2 ) 

Z R (q 2 ^ 2 ) q 2 + M 2 (q 2 ) 



2(p 2 q 2 ~(p-q) 2 ) 



u° M ( P 2 ,i 2 y), 



(29) 
(30) 

(31) 
(32) 



Let us verify that this truncation preserves the multiplicative renormalizability of the solutions. From Eq. (|T(j) one 
can show that multiplicative renormalization is satisfied if solutions renormalized at /i 2 and v 2 arc related by: 



Z R (P 2 ^ 2 ) 



Zr( P 2 ,H 2 ) 
Zr{v 2 ^ 2 ) 



Z 2 (v 2 ,A 2 ) = Z R (v 2 ,v 2 )Z 2 (p 2 ,A 2 ). 



(33) 



Now assume one has found solutions Z R (p 2 ,^ 2 ), M(p 2 ) of Eqs. (p9|), (j30l), for all p 2 , renormalized at a scale fi 2 , i.e. 
with Z 2 (^ 2 , A 2 ) such that Zr(^ 2 ,// 2 ) = 1. If one multiplies both equations (p9|), ( |30| ) by Zr(v 2 , /i 2 ), it is clear that 
Zr(p 2 , v 2 ) of Eq. ([33j), and the original mass function M(p 2 ) are solutions of the new equations, satisfying Zr(v 2 , v 2 ) — 
1, and with renormalization constant Z 2 (v 2 , A 2 ) obeying Eq. (63). The ability to construct the solutions renormalized 
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at v 2 starting from the solutions renormalized at /i 2 , following the principles of multiplicative renormalization defined 
in Eq. (|Io|), demonstrates that the novel truncation respects the multiplicative renormalizability of the solutions. Such 
a construction is not possible for Eqs. (0), ( |l8| ) of the bare Abelian approximation, hence showing its violation of 
multiplicative renormalizability. 

Also note that the solutions of the truncated equations (p9|), ( |30| ) satisfy the leading order resummed perturbative 
results. This feature was demonstrated earlier for the Abelian case jy], ^9), and even though the newly proposed 
equations differ from the former in their Z-dependence, their UV-limits are identical, as Zr(p 2 ) = 1 to leading order 
in the Landau gauge, and hence, their mass functions have the same UV behavior. For systems where Zji(p 2 ) ^ 1 the 
Abelian approximation does not reproduce the correct perturbative behavior for the quark functions, while truncations 
similar to Eqs. (|27j), (^8|) still do. This has previously been shown for the gluon and ghost dressing functions in QCD 
p5|, which obey similar equations. Also note that a truncation of the fermion equation in QED 3 |30[ ], constructed in 
order to reproduce the correct perturbative behavior, features a similar linearity in the dressing function Z which, 
as was shown above, is sufficient to achieve multiplicative renormalizability. The treatment of QCD is, however, 
significantly different from QED, as the runnin g co upling only enters the equations after introducing appropriate 
ghost field corrections as proposed in Eqs. (P7j), (P8|). 

The fact that the simple truncation ( p7| ) , (Eq) satisfies these two important features of gauge theories is not a mere 
accident. Similarly to our work on the gluon equation in Ref. p5f , the cancellations assumed in our Ansatz can be 
related to the cancellation of the quantum corrections coming from both the full quark-gluon vertex and the integral 
over the DSE kernel, in a way similar to Mandelstam's work on the gluon equation fll4f . Indeed, the quark-gluon 
vertex Slavnov- Taylor identity, which is a consequence of gauge invariance, is 

{p-krY^k ] p-k)^G{(p^k) 2 ){S F 1 {k)-S F \p)) , (34) 

which shows that the bare vertex receives nonperturbative corrections proportional to G/Z. Although the Slavnov- 
Taylor identity does not constrain the transverse part of the vertex, multiplicative renormalizability requires the 
complete full vertex to have corrections of a similar nature. Furthermore, Mandelstam also showed in Ref. (T^j how 
perturbative loop corrections to the propagators in fact introduce a double G/Z correction, and we herein assume 
that this is not just a perturbative feature, but rather a true nonperturbative one, engendered by the integrations 
over the full Green's functions in the self-energy integrals. Hence, in our analysis, one factor Z/G is assumed to 
cancel the quantum corrections of the full vertex, and the other factor Z/ G to cancel the nonperturbative corrections 
generated by the integrals over the kernels. This nonperturbative Ansatz is consistent with perturbation theory at 
large momenta, hence yielding the correct resummed perturbative behavior of the quark functions, but moreover it 
preserves the multiplicative renormalizability of the solutions for any choice of renormalization point. Note that, even 
though the full vertex has not been constructed explicitly, it is genuinely nonperturbative in our truncation - however, 
the final equations contain bare vertices only, as a consequence of the non-trivial cancellations assumed in our Ansatz. 

The main aim of this study is to investigate the realization of dynamical chiral symmetry breaking in the proposed 
MR preserving truncation, where the kernel is explicitly driven by the running coupling. From the derivation of Eqs. 
(p9|), j30| ) we note that the running coupling in these equations is the true QCD running coupling, and not an effective 
coupling as was the case in the Abelian approximation, Eqs. (|l7j), (|l8|). The running coupling connects the quark 
equation to the gauge sector, and in the next section we discuss the strong running coupling in more detail. 



III. THE STRONG RUNNING COUPLING 

Formally the truncated quark equations ( p9f ) , ( |30| ) have to be solved in a coupled system together with the equations 
for the gluon and ghost propagators presented in Ref. J|^] , such that the back reaction of the quarks on the gluon 
vacuum polarization be taken into account correctly in a truncation that preserves multiplicative renormalizability. 
However, as shown in Ref. p5| , the pure gauge gluon-ghost system itself is not yet fully understood, and we can 
circumvent this complicated analysis, and postpone the study of the coupled equations to a later stage, by introducing 
a realistic model for the running coupling, thus decoupling the quark equations from the gluon-ghost system. 

We will construct a model running coupling using information originating from perturbation theory, from Dyson- 
Schwinger studies of the coupled ghost-gluon system, and from recent lattice calculations. 
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Asymptotic behavior 



Renormalization group equation improved leading order perturbation theory yields the asymptotic ultraviolet be- 
havior of the running coupling, including the quark contribution to the vacuum polarization: 



where 



«(9 2 ) 



-in 



f3 \og{q 2 /K 2 



qcd; 



HN r - 2N 



f 



(35) 



(36) 



with N c colors and Nf flavors, and where the leading order Aqcd can be determined from some high energy experiment. 
Note that this leading order resummed behavior is also reproduced in the solutions of the nonperturbative DSEs, when 
computed using the truncation described in Ref. [25j| , which preserves multiplicative renormalizability. 

From an analysis of the coupled ghost and gluon DSE system, von Smekal et al. |18 showed that, in a specific 
truncation scheme, the gluon and ghost propagators obey power laws in the infrared region, and the running coupling 
has an infrared fixed point. Furthermore, the leading infrared behavior of the gluon equation is completely determined 
by the ghost loop contribution to the vacuum polarization. These results were confirmed by studies of Atkinson and 
Bloch [2(| using different vertex Ansatze and angular approximations. However, the various truncations used in 
these studies all violate the principle of multiplicative renormalizability, and Ref. p?| proposed a truncation where 
multiplicative renormalizability is respected, and which possesses the distinctive feature that all diagrams contribute 
to the leading infrared behavior of the gluon equation. A recent investigation of the contributions of the two-loop 
diagrams to the gluon vacuum polarization (26) in this MR truncation reveals a relatively large infrared contribution of 
the four-gluon or squint diagram which is essential, when combined with the ghost and gluon loop contributions [ p5[ , 
to ensure the existence of the propagator power laws, and of an infrared fixed point for the strong running coupling. 
This is a satisfying observation as the infrared power law solutions for both the gluon and the ghost propagators, and 
the infrared fixed point of the strong coupling are also supported by a number of recent lattice calculations pil B2, E3[ . 



The model 



Taking into account all this evidence, we herein build a model running coupling which has an infrared fixed point, 
and a constraint on its value can be d erive d from the ghost DSE only, as shown in Appendix [X], independently of the 
details of the gluon DSE. From Eq. (All) we compute that, for SU(3), the fixed poin t satisfies 2ir/3 < «o < 47r/3, 
depending on the precise value of the exponent n of the propagator power laws (A3), and combining with recent 
lattice results giving k ~ 0.5 [Ell E2], ^3), we assume a preferred value ao = 5ir/6 ~ 2.6. 

Even though we pinned down the leading asymptotic behavior of the running coupling at large and small momenta, 
we still need to know the behavior of the running coupling for all momenta in order to solve the quark equations (|2S|), 
(|30|). To model the intermediate momentum region, where the coupling evolves from the infrared fixed point into 
the logarithmic tail, we use information collected from the study, Ref. pj| , of the coupled ghost-gluon system. The 
numerical results and the analytically calculated infrared asymptotic series in that study showed that the coupling 
remains very close to the infrared fixed point almost up to q 2 = Aq CD , crosses the leading order perturbative curve, 
and thereafter drops quite rapidly to rejoin the perturbative curve, from above, around q 2 = 10 Aq CD . 

We propose a functional expression for the running coupling which satisfies the various constraints mentioned 
above, making use of rational polynomial approximations multiplying the required asymptotic behaviors. In its 
simplest form our approximation contains only one, dimensionless, parameter cq apart from the infrared fixed point 
«o and the intrinsic scale Aqcd, 



a(q 2 ) = a(t A 2 



QCDJ 



C{) 



f 2 



c a 



1 



1 



4-7T 



A, \\ogt t-i 



(37) 



where t = <7 2 /Aq CD j|3). The first term in the square brackets of Eq. ( p7\ ) is responsible for the infrared fixed point, 
while the term between round brackets yields the correct leading order logarithmic tail, where the simple pole at 
Aqcd has been subtracted to make the coupling analytic for all spacelike momenta [j3l| . To match the intermediate 
region with the results of Ref. jl9j, we choose Co = 15, and we illustrate this model for the running coupling in Fig. 
[I]. Note that in contrast to the effective coupling often used in the application of DSE to hadron phenomenology, the 
running coupling (j37|) does not have any infrared enhancement, and can be considered as a smoothened version of the 
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10000 



FIG. 1: Running coupling a(q 2 ), as given by Eq. ( J37| ) with ao = 2.6 and Co = 15, and q 2 given in units Aq Cd . The coupling 
has an infrared fixed point ao, and an ultraviolet behavior consistent with resummed leading order perturbation theory. 



concatenation of the infrared fixed point with the leading order logarithmic tail, where the smoothcning is performed 
in a way suggested by previous studies of the gauge sector. The meaning and determination of the intrinsic QCD 
scale Aqcd in the model ([57]), will be discussed in the result section. 



IV. NUMERICAL RESULTS 



We now study how dynamical chiral symmetry breaking gets realized in Eqs. (|29|), ( p0| ) with kernels driven by 
a(q 2 ) of Eq. (37). Note that a(q 2 ) possesses an intrinsic scale Aqcd, with which all momenta are scaled, and hence 
all momenta and masses in the quark equations can be expressed hereinafter in units of Aqcd- The equations are 
solved numerically for the chiral case, and the results, which satisfy MR, are compared with those computed in the 
MR-violating Abelian approximation. We investigate the sensitivity of the generated mass scale to the parameters 
of the model coupling, and discuss the determination of Aqcd of the model. We also show the results of the quark 
equation in the massive case. 

The seeds of Eq. (|9|), (^) depend on the unknown renormalization constant Z-i and bare mass mo, and they 
are eliminated in the usual way by subtracting each equation at the renormalization point /z 2 from that at p 2 , and 
imposing the renormalization conditions Z R (fi 2 ) — 1, M(fi 2 ) = m^. All integrals in the subtracted equations are 
finite, and the ultraviolet regulator can be taken to infinity. After introducing spherical coordinates, and performing 
two trivial angular integrations, we find 



1 



Zr{x) 
M(x) 



1 



Cf_ 

^ 2 J 

3C F 



dy 



V 



Z R (y)y + M 2 (y) J 
M(y) y 



d9 sin 2 e 



a(z) 



3y/yxcos9 2xysiv?9 1 



dy 



d9 sin 2 e 



a(z) 



— (x <-> p 2 ) 



,(38) 
(39) 



Z R {x) 11 2ir 2 J y Z R {y)y + M 2 {y) J 

with x = p 2 , y = q 2 and z = x + y — 2^fxycos 8. Note that from now on the renormalization point dependence of Zr 
will be implicitly understood, Zr{x) = Zr(x,[i 2 ). 



Chiral case 

In the chiral case the mass equation ( |30| ) is homogeneous, as ^"^o = 0, and no subtraction is needed for the mass 
equation - Eq. (|3^) is replaced by 

Mix) _sc F r dy mp_ V r de sin2 e «w . (40) 



Z R (x) 2tt 2 J Z R (y)y + M 2 (y) J 



We substitute our model running coupling (37) in the quark equations (p8|), (ktfj) and compute the solutions of 



these equations using a numerical procedure first constructed to study dynamical chiral symmetry breaking in strong 
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FIG. 2: Solutions for the mass function M(x) and quark dressing function Zr(x), from Eqs. (^), (0), renormalized at 
fi 2 = Aq CD , using a(x) given by Eq. (pS7|), with ao = 2.6 and co = 15, in the chiral case. Preservation of multiplicative 
renormalizability means that a different choice of renormalization point will merely multiply Zr{x) by a constant factor, and 
leave M(x) unchanged. 



coupling QED jj2| , and later also applied to solve the coupled gluon-ghost equations |T9| . The method uses Chebyshev 
polynomials to approximate the unknown functions, and solves for the polynomial coefficients using the quadratically 
convergent Newton iterative method. 

The evolution of M(x) and Zr(x) as a function of momentum is shown in Fig. |^, for 10~ 5 < x/Aq CD < 10 5 , and 
renormalization point fi 2 = Aq CD . The mass function is expressed in units of Aqcd, which sets the scale in the running 
coupling (|37]). Note that, in our MR preserving truncation, a different choice of renormalization point will merely 
multiply the whole function Zr(x) by a finite, constant factor over the complete momentum range, and will leave the 
mass function M(x) unchanged. The strength of the dynamical chiral symmetry breaking can be characterized by 
the value of the mass function at the origin: 

M(0)» 1.057 Aqcd- (41) 

The generated quark mass is of the order of the extension of the infrared plateau of the coupling, and the absence of 
fine-tuning in the coupling, to obtain such a mass scale, points to a certain naturalness in the MR truncation. 

In the chiral case the ultraviolet behavior of the dynamical mass function satisfies the well-known asymptotic 
formula H,||: 

M(x) ~ 2 ^ lm ~ {qq) ° (42) 

with 7 m = 12/(33— 2Nf). A large momentum fit yields the value for the renormalization-point independent vacuum 
quark condensate, which is also representative for the strength of dynamical chiral symmetry breaking: 

- (qg)° « (0.70 Aqcd) 3 - (43) 



Abelian approximation 



An important part of our analysis is to evaluate the influence of the multiplicatively renormalizable approach on 
the dynamically generated mass scale. The only difference between the MR Eqs. (^9|), ( j30|) and Eqs. (17), (lj) of the 
bare vertex, Abelian approximation, is the factor Zn(y) which appears in the denominator of the integration kernels 
in the former, while it is found in their numerator in the latter, in a way violating MR. Using our model running 
coupling (37) as effective coupling in the Abelian approximation, we find a dynamically generated mass that strongly 
depends on the renormalization point, as can be seen in Table |[ The value which is usually quoted in Abelian studies 
is that of the unrenormalized case, for which the dynamical mass is M(0) ~ 0.326 Aqcd, about three times smaller 
than in the herein presented MR preserving truncation. Our MR truncation clearly disagrees with the use of a bare 
quark-gluon vertex, and the implicit dressing of the vertex assumed in the MR truncation is essential. The violation 
of MR, and hence the renormalization point dependence of the generated mass scale makes the bare vertex Abelian 
truncation ill defined. 
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2 

M 


M(O) 


MR 


any 


1.057 


Abelian 


unren 


0.326 




10 5 


0.327 




10 3 


0.355 




1 


1.262 




1(T 3 


1.398 




10~ 5 


1.398 



TABLE I: Dynamically generated mass M(0), in 
units Aqcd, versus renormalization point (in 
Aq CD ) in MR- violating Abelian approximation, 
compared with MR solution. 



A?r„ 

il UV 


M(0) 


0.1 


0.076 


1 


0.562 


10 


1.044 


100 


1.057 


10 3 


1.057 



TABLE II: Dynamically generated mass M(0), in 
units Aqcd, versus upper integration limit A.fj V 
(in A QCD ) in Eq. ©. 



We furthermore observe that the mass scale generated in the MR truncation depends very little on the momentum 
evolution of Z. To show this, it suffices to decouple the mass equation ( pio| ) from the Z-equation ( p8[ ) by forcing 
Zr(x) = 1 for all x in the mass equation. The generated mass is then M(0) = 1.082 Aqcd, only 2% off the value of 
the coupled system of equations. 



Model sensitivity 



We now study the sensitivity of the solutions of the MR truncated quark equations to the features of the model 
( |37| ) for the running coupling. 

First, we examine which part of the kinematical regime of the running coupling (37) mainly contributes to the 
dynamical mass in the MR truncation. We substitute the solutions to Eqs. (p8|), (^0|), shown in Fig. |^, in the right 
hand side of Eq. (M), and calculate the mass integral for varying values of the upper integration limit Ay V . The 
results are given in Table ||, and we find that slightly more than half of the dynamical mass is generated by the 
infrared fixed point region of the coupling, p 2 < Aq CD . Most of the remaining mass is generated in the transition 
region, Aq CD < p 2 < 10Aq CD , while the logarithmic tail only contributes about 1%. Hence, the dynamical mass is 
mainly generated in the nonperturbative region of the coupling, and is virtually independent of its precise ultraviolet 
behavior. It is therefore not essential, in the context of dynamical chiral symmetry breaking, that the ultraviolet 
behavior of the coupling ( |37| ) only reproduces the RGE improved leading order perturbative results. 

Next, we investigate the sensitivity of the results to the parameters of our model coupling (|37|). As shown in 
Appendix the value of the fixed point depends on the details of the infrared analysis of the gluon-ghost DSE 
system, with bounds 27r/3 < ao < 47r/3, and preferred value 5tt/6. In Fig. [| we show how the generated mass changes 
as we vary the value of the fixed point over a slightly wider range, 1 < ao < 12, and some of the main values are 
also tabulated in Table [II. In the figure we also plot the corresponding results for the unrenormalized, bare vertex, 
Abelian case discussed earlier. As could be expected there is a critical value a™ 1 ' « 0.90 (for Co = 15) below which 
chiral symmetry is not spontaneously broken. However, the fixed point value dictated by the QCD equations is well 
above this value, and the quarks do acquire an appropriately sized dynamical mass in the MR truncation. 

When constructing the model coupling (|37]), the shape of its drop in the intermediate region, Aq CD < p 2 < 10 Aq CD , 
where the coupling flows from the infrared asymptotic behavior into its perturbative logarithmic tail, was inferred 



ao 


M(0) 


1.0 


0.028 


1.8 


0.548 


2.6 


1.057 


3.4 


1.494 


4.2 


1.874 



TABLE III: Sensitivity of the dynamically gen- 
erated mass M(0), in units Aqcd, to changes in 
the infrared fixed point of the coupling, ao of Eq. 
@, with c = 15. 



CO 


M(0) 


5 


0.824 


10 


0.964 


15 


1.057 


20 


1.129 


40 


1.325 



TABLE IV: Sensitivity of the dynamically gener- 
ated mass M(0), in units Aqcd, to changes in the 
parameter Co of Eq. with ao = 2.6. 
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FIG. 3: Variation of the generated mass M(0), in units of 
Aqcd, with varying cto in Eq. (jj^), and Co = 15, for our 
MR truncation scheme, and in the Abelian approxima- 
tion. The vertical line shows the preferred value ao = 2.6. 
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x [ a qcd1 



FIG. 4: Running coupling a(x), given by Eq. (^) with 
ao = 2.6, for a range of co's. To focus on the influence 
of Co on the intermediate momentum region, the plotted 
momentum range is taken 0.1 < £/Aq CD < 1000. 



from the numerical results of Ref. [jl9| . Even though we are confident that the qualitative features of this study will 
remain valid outside the scope of its specific truncations (angular approximation, violation of MR, disregard of quark 
loop and two- loop diagrams), more will be learned about the intermediate regime when solving the coupled quark- 
gluon-ghost set of equations in an MR preserving truncation, as we plan to do in future work. However, the sensitivity 
of the generated mass scale with respect to a change in intermediate behavior can already be studied by varying the 
parameter cq of our model coupling (|37|). The behavior of the running coupling, as cq is varied, is illustrated in Fig. ||, 
and the dynamically generated mass is tabulated in Table IV, for selected values of Cq. The parameter c is related to 
the position of the transition region, where the running coupling drops from its fixed point to rejoin the perturbative 
curve, and the dynamical mass gradually increases with increasing cq. 

Note that another approach yielding a finite infrared value of the running coupling can be derived from an analytic 
continuation of the perturbative coupling as shown in Ref. |3lJ]. In contrast to the model used above, the coupling 
derived in that approach has no infrared plateau, and is exactly given by expression ( |37| ) with cq = 0, yielding a finite 
infrared value a(0) = 1.51, for Nf — 4. This running coupling does not yield enough integration strength in the gap 
equation, as the computed dynamical fermion mass M(0) = 0.036 Aqcd is clearly too small for hadronic applications. 



Determining Aqcd 

As the running coupling a{q 2 ) of Eq. ( ]37| ) is a function of g 2 /Aq CD , all mass scales in the quark equations are 
expressed in units Aqcd, and the value of Aqcd has to be determined by matching theory with experiment. 

Determining Aqcd of our model by matching a in the perturbative region is not a reliable procedure, as the 
perturbative behavior of our model coupling ([37]) only agrees with perturbation theory up to resummed leading order, 
and small changes in a in the logarithmic region will have a big impact on the value of Aqcd, and on all related 
hadronic masses, as can be inferred from Eqs. ([l"l|), ([f3"|). 

In perturbation theory Aqcd corresponds to the position of the pole in the logarithmic expansion of the running 
coupling, and depends on the renormalization scheme, and on the mass threshold matching procedure p5[ . In the MS 
scheme its five flavor value determined from experiment is A^jL = 208^23 MeV [^H). However, Aqcd changes as flavor 
thresholds are crossed, and the running of the perturbative coupling should be calculated using the appropriate Aqq^. 
Taking into account mass thresholds, one finds A— = 309 + qn MeV below the b-quark mass, and A— = 377 + ™ MeV 
below the c-quark threshold (as long as perturbation theory remains valid). 

The nonperturbative model ( |37| ) has no perturbative pole, and there Aqcd is uniquely related to the evolution of 
the coupling (^) , and corresponds to the scale where the coupling starts its descent off the infrared plateau towards 
the perturbative regime (see Fig. [I]). As this definition of Aqcd is at variance with the perturbative definition, a 
straightforward identification of its value with the perturbative values is not possible. Our aim herein is to investigate 
nonperturbative effects, like dynamical chiral symmetry breaking, hence Aqcd of our model should be determined by 
matching some nonperturbative quantity related to the scale of hadronic physics. 
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FIG. 5: Comparison of the Dyson-Schwinger solutions for the mass function M(x) and quark dressing function Zr(x) in the 
chiral case, with the lattice results of Ref. |37| (chiral extrapolation). The Dyson-Schwinger solutions are computed from Eqs. 
©, @, using a(x) gwen by Eq. @, with a = 2.6, c = 15, and with Aqcd = 330 MeV and 550 MeV (see text). The 
lattice results of Ref. [|37[ are plotted as a function of the kinematic lattice momentum. 



If we take Aqcd = 330 MeV in the model coupling (p7|), the MR truncation generates a dynamical quark mass 
M(0) « 350 MeV and a chiral condensate -(qq}° ~ (230 MeV) 3 , according to Eqs. ([H]), (g), which are of the right 
order of magnitude for hadron phenomenology. Moreover, this value of Aqcd is consistent with the above mentioned 
perturbative values, taking into account their different definitions, as the coupling grows strong in the same momentum 
window. 

Furthermore, for such a value of Aqcd, our model coupling (|37j), whose ultraviolet behavior is only modelled up to 
leading order in perturbation theory, yields a perturbative value a(Mz) = 0.134 (Mz = 91.187 GeV). Particle physics 
phenomenology gives a world-average a s {Mz) = 0.118 for the strong coupling |3j|, and the leading order contribution 
to this value, computed from perturbation theory, is al~ loop (Mz) = 0.135, including mass-threshold effects Q. We 
can thus conclude that our model coupling produces results which are consistent with both nonperturbative and 
perturbative physics. 

The numerical results obtained herein can also be compared with some of the latest lattice results on the quark 
propagator in the Landau gauge (with Nf — 0) pTj]. Fig. || shows the lattice data, and the DSE results with 
the phcnomcnologically preferred value Aqcd = 330 MeV, and also with Aqcd = 550 MeV, which yields a better 
agreement with the lattice points in the intermediate momentum regime (3 GeV 2 < x < 50 GeV 2 ). Note that below 
3 GeV 2 the lattice results have a slower increase of the mass function, and a steeper decrease of the dressing function. 
It is this discrepancy at small momenta that affects the value of Aqcd when matching the DSE-results with the lattice 
data. Possible origins are lattice artefacts, vertex truncations, or inaccuracies in the modeling of the coupling (p7|) in 
the intermediate momentum regime. As the size of the dynamically generated mass is crucial for hadronic physics, 
this disparity should be investigated further. 



Massive case 



Next we extend our analysis to the massive case, ^"t-o 7^ in the mass equation (p0| 



in Eq. (p9[), where 



n chiral 



or in practice > m^ hlral 

= M ciural (/i 2 ). In Fig. |6j we plot the mass and quark dressing functions for varying 
renormalized mass, to^/Aqcd = 0.0001,0.001,0.01,0.1,1, together with the chiral case (with renormalization point 
// 2 /Aq CD = 10 5 ). In the massive case the logarithmic behavior satisfies the RGE improved perturbative result, 



M(p 2 ) = m p 



a(p 2 ) 



log(pVA 2 QCD ) 



(44) 



where rh = m M [log(/Lt 2 /Aq CD )] 7 " 1 , and 7 m = 12/(33 — 2Nf), which is the correct anomalous dimension of the mass 
function. We observe that the log-tail of the mass function sets in for ever smaller momenta as the renormalized mass 
is tuned above the chiral value. Note again that the preservation of multiplicative renormalizability allows us to choose 
an arbitrary renormalization point, without altering the physical content of the results: changing the renormalization 
point merely multiplies Zr(x) by a constant factor at all momenta, and leaves the mass function M{x) unchanged. 
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X [Aqcd] x[a£ cd ] 

FIG. 6: Mass function and quark dressing function solutions from Eqs. (^), (^), using a(x) given by Eq. (^), with ao = 2.6 
and Co = 15, for various values of renormalized quark mass m Ml with renormalization point fi 2 /Aq CD = 10 5 . Note the influence 
of the quark mass on the onset of the logarithmic tail in M(x). 

Discussion 

As shown before, the Abelian approximation, which violates MR, requires a strong infrared enhancement of the 
effective coupling to achieve a strong enough breaking of the chiral symmetry, and such an enhancement seems 
inconsistent with the characteristics of the gauge sector. A first attempt to extend the DSE studies beyond the 
Abelian approximation by solving a truncated set of DSEs for the quark, gluon and ghost propagators simultaneously 
p4j finds a dynamically generated chiral condensate of (122 MeV) 3 (with Nf = 4), which is about a factor two too 
small for phenomenological purposes, apparently confirming the problem to find the necessary integration strength 
in the quark equation. The study uses a QCD analog Q of the quenched QED Curtis-Pennington vertex |59| in an 
attempt to preserve gauge invariance and MR. This is however not sufficient: the vertex renormalization constant 
Zif is constructed such that the solutions to the equations have the correct leading order perturbative behavior, but 
it is not consistent with the chosen nonperturbative vertex Ansatz, and its renormalization point dependence violates 
multiplicative renormalizability. 

In Rcf. a method was presented to solve the renormalized quark equations in an MR truncation of quenched 
QED. There the explicit construction of Z\f is circumvented by eliminating the renormalization constants from the 
equations. It is however not clear if and how the method could be extended to QCD. Moreover, in contrast to our 
MR truncation, their method requires an explicit construction of a nonperturbative vertex Ansatz satisfying the 
Ward-Takahashi identity and conditions of MR, and the two truncation schemes are not equivalent. 

The MR truncation scheme presented in this paper leads to solutions that satisfy the principles of multiplicative 
renormalizability in an elegant way, without explicit constructions of full vertices or renormalization constants, while 
naturally producing a large enough dynamical breaking of chiral symmetry in QCD. In a next step, the MR truncation 
derived for the gluon-ghost system pf| should be combined with that proposed herein, to solve the coupled set of 
DSEs for quark, gluon and ghost propagators, thus avoiding the explicit construction of a model running coupling. 
Consequently, an MR preserving truncation scheme for the Bethe-Salpeter equation, describing hadronic bound states, 
should be constructed in a way consistent with the quark equation derived herein, which respects the Goldstone boson 
nature of the pion. 

V. CONCLUSIONS 

We have reformulated the coupled set of continuum equations for the quark dressing function and dynamical mass 
function in QCD, such that, in the Landau gauge, all renormalization constants are eliminated, and the multiplica- 
tive renormalizability of the quark dressing function, and renormalization point invariance of the mass function are 
manifest. The kernels of the equations explicitly depend on the momentum evolution of the strong running coupling. 
The new formulation allowed us to construct a truncation which preserves these properties, and also satisfies the 
leading order renormalization group equation improved perturbative results. This is achieved by assuming non-trivial 
cancellations involving the full quark-gluon vertex in the self-energy loop, such that no explicit construction of the 
full vertex and vertex renormalization constant are needed. 
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In order to perform numerical calculations of the dynamically generated mass scale, without solving the coupled 
quark-gluon-ghost DSE system, we decoupled the quark equations from the gluon and ghost equations by introducing 
a model running coupling. This coupling satisfies the leading order ultraviolet behavior known from perturbative 
QCD, and its nonperturbative behavior was constrained using the most recent results of gluon-ghost DSE and lattice 
studies. The coupling has an infrared fixed point («o ~ 2.6), and its behavior in the intermediate momentum region 
reproduces the qualitative features of the numerical DSE studies. 

The integral equations for the quark propagator are solved numerically for the chiral and massive cases, and the 
fermion mass and chiral condensate generated dynamically in the MR preserving truncation is of the order of the 
extension of the infrared plateau of the coupling. The generated scale is about three times larger than what is found 
using the MR-violating, bare vertex, Abelian approximation, and of the right order of magnitude to perform hadron 
phenomenology, without infrared enhancement in the strong coupling. 
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APPENDIX A: INFRARED BEHAVIOR OF THE STRONG RUNNING COUPLING 



We briefly show how the value of the infrared fixed point ao of the strong coupling should be determined from the 
coupled Dyson-Schwinger equations for the gluon and ghost propagators. The contributions of the various diagrams 
in the gluon equation are however not yet fully understood, and we describe how the ghost equation alone can give 
us tight bounds on the value of ao- 

The DSE for the gluon propagator is 

[zvwr 1 = [dim 1 - - - 7$ (?o - 4up) - - km , (ai) 

where D^ v and are the bare and full gluon propagators, and the vacuum polarization includes contributions from 
the ghost loop, gluon loop, three-gluon diagram, four-gluon diagram, tadpole diagram and quark loop. The DSE for 
the ghost propagator, in Euclidean space, is 

[A(p)]- X = [A°(p)] _1 - N c g 2 J 0^Gl(p,q)A( q )G u (q,p)D^(r), (A2) 

where N c is the number of colors, go is the bare coupling, A and A are the bare and full ghost propagators, G° 
and Gfj, the bare and full ghost-gluon vertices, and r = p — q. The general expressions for the full gluon and ghost 
propagators, defining the gluon and ghost dressing functions F(p 2 ) and G(p 2 \ were given in Eqs. (Q), (||). 

As mentioned in the main body of this paper, the coupled equations Eqs. (Al), ( A2) have been solved simultaneously 
in a number of previous studies using various truncations, and there is convincing evidence that F and G obey power 
laws in the infrared, when all the diagrams in the gluon equation are taken into account properly: 

F R (x, fj?) ~ x 2K , G R {x^ 2 )^x- K , (A3) 

with x = p 2 , and where Fr and Gr are the renormalized gluon and ghost dressing functions defined in Eq. (|l3|). 
Because of Eq. (|23), these power laws lead to an infrared fixed point for the running coupling, 

a a = lim a(fi 2 ) Fr(x, /j 2 ) G r (x, ^l 2 ) — ► constant . (A4) 

x— >0 

As explained in detail in Ref . pq] , the ghost and gluon equations each typically yield a relation between the infrared 
fixed point ao and the leading infrared exponent k: 

Nc a = — ' N " a o = — ' ( A5 ) 
Xgh[K) Xgl{K) 



when substituting the power laws (A3) in the coupled equations ( |Al[ ), (|A2j) , and equating the coefficients of the 
leading power of x for x — > on both sides of the equations. The functions Xgh{ K )i Xgi( K ) are computed from the 
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vacu um polarization integrals, and a consistent infrared power solution requires the gluon and ghost expressions in 
Eq. (|A5| ) to be satisfied simultaneously: 



Xgh{n) = Xglfa)- 



(A6) 



The solution of this equation yields the value of the leading infrared exponent re, and the corresponding ao can then 
be computed from E q. (Aq ). 

In a recent study p5pwe have shown that the gluon function Xgi( K ) is a complicated object as it requires the 
computation of the nonpcrturbative loop integrals for all the gluon vacuum polarization diagrams in the infrared. 
However, if the power laws are valid, then a consistent value of re exists, and the ghost equation can, on its own, 
provide very useful information. Indeed, the first identity of Eq. ( A5) rel ates the value of the infrared fixed point ao of 



the strong coupling to the infrared exponent re of the dressing functions (A3), without needing additional i nfor mation 
from the gluon identity. Therefore, we now further analyze the infrared behavior of the ghost equation (A2) in the 
bare ghost-gluon vertex approximation, and compute Xgh ■ 



In the Landau gauge, the bare vertex approximation of Eq. (A2) for the renormalized ghost dressing function 
becomes, 



G R (x,fi 2 ) 



d A q T (x, y, z) G R {y, fi 2 ) F R (z, fi 2 ) , 



(A7) 



where we introduced the full propagators (|4|), (^), the renormalized quantities (p"2|), (p"3|), used Z\ = \ (in the Landau 



gauge), set x = p z , y 



2 

q , z 



r 2 , and 



7b( w ) = -|--2 + ! — 2 



1 " 

x 2z 



Axy 



(A8) 



After substitution of the power laws (A3) in Eq. (A7), the right hand side yields a sum of integrals of the form 

d A q x a / z~< , (A9) 



with a + (3 + j = K — 2. In Ref. j2Q] these integrals were solved using spherical coordinates, and the results were 
expressed in terms of generalized hypergeometric functions. A more concise, equivalent expression can be derived 
plj by noting that integrals of type (A9) are typically encountered when applying dimensional regularization in 
pcrturbative calculations, and they are readily computed by introducing Feynman parameters, yielding |43] 



I(a,b) = 



1 



d A q 
(2tt) 4 y a z b 



1 



16tt 2 



r(2-a)r(2-6)r(a + 6-2) a;2 _ a _ fe 



r(a)r(6)r(4 - a - b) 



Both sides of Eq. (A7) yield a leading infrared power x K , and equating their coefficients gives 

2tt r(3 - 2re)r(3 + re)r(l + re) 



3N C 



T 2 (2- re)r(2re) 



(AlO) 



(AH) 



where the Ga mma function recurrence relation was used repeatedly to bring the expression in its simplest form. 

Expression (All) gives the relation between the infrared fixed point of the coupling and the infrared exponent re 
imposed by the ghost equation, and we illustrate this dependence for SU(3) in Fig. ^. As mentioned before, the gluon 
equation will be consistent with this expression for some value of re if the infrared power laws for the propagators are 
valid. Previous DSE studi es |18|, jl9| , |20| , ^l[ yielded values for re between 0.4 and 1.0, depending on the truncation, and 
recent lattice calculations 21 , |22]| seem to indicate that re is close to 0.5. Although the exact value of re depends on the 
consistency condition ( |A6|) between the ghost and gluon DSEs, we see from Fig. that the value of ao is tightly bound, 
and will lie in the interval [27r/3, 4tt/3] for SU(3), and the preferred value used herein is o.o(k — 0.5) = 57r/6 « 2.6. 
Note that the large values for the fixed point, ao ~ 10 — 12, mentioned in some previous studies JlJ, have been 
shown to be artefacts of angular approximations in the integration kernels pfj . 
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